Off Equilibrium Study of the Fluctuation-Dissipation Relation in the Easy- Axis 
Heisenberg Antiferromagnet on the Kagome Lattice 
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Violation of the fluctuation-dissipation theorem (FDT) in a frustrated Heisenberg model on the 
Kagome lattice is investigated using Monte Carlo simulations. The model exhibits glassy behaviour 
at low temperatures accompanied by very slow dynamics. Both the spin-spin autocorrelation func- 
tion and the response to an external magnetic fleld are studied. Clear evidence of a constant value 
of the fluctuation dissipation ratio and long range memory effects are observed for the flrst time in 
this model. The breakdown of the FDT in the glassy phase follows the predictions of the mean field 
theory for spin glasses with one-step replica symmetry breaking. 

PACS numbers: 64.70.Pf, 75.40.Gb, 75.40.Mg 
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Glassy systems are characterized by long relaxation 
times which increase at low temperatures and eventu- 
ally exceed experimental observation times. At higher 
temperatures, the systems are ergodic and true equi- 
librium can be reached. However, below the glass 
transition temperature, Tg, the systems can become 
trapped in metastable states and often display strong 
non-equilibrium effects such as aging. The term aging 
refers to the fact that both static and dynamic properties 
depend on the time interval after which the system has 
been quenched into some non-equilibrium state. Struc- 
tural glasses Q and spin glasses |3| exhibit these features 
and are essentially out of equilibrium on experimental 
time scales. In the case of spin systems, this behaviour 
is usually found in sytems with some type of quenched 
disorder. The disorder often imposes severe restrictions 
on the spin rearrangements and the evolution towards 
equilibrium requires many degrees of freedom to act co- 
operatively. The complexity of the phase space can lead 
to the violation of the fluctuation dissipation theorem 
(FDT) below a characteristic temperature. 

In the last decade off equilibrium approaches have been 
used quite succesfuUyto describe systems which show a 
very slow dynamics [2|. A typical example is the zero- 
field cooling ^ experiment, in which a sample is cooled 
in zero magnetic field to a low temperature at time t — 0. 
After a waiting time a small magnetic field is applied 
and subsequently the time evolution of the magnetiza- 
tion is recorded. It is often observed that the relaxation 
of the magnetization becomes slower as the waiting time 
tu, is increased. The system never reaches thermody- 
namic equilibrium and, consequently, time translational 
invariance (TTI) and the fluctuation dissipation theorem 
are no longer valid. However, the system is in a quasi- 
equilibrium state for time scales which are much longer 
than microscopic time scales but shorter than the relax- 
tional time scale of the slowly relaxing degrees of freedom. 



In this glassy regime a generalized form of the FDT has 
been proposed 
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where C (t + , tw) is a double time correlation function 
and R{t+tw, t^) with i > is the associated conjugate re- 
sponse function. In this relation, T is the heat bath tem- 
perature and the function X {t + tw , tw) is called the fluc- 
tuation dissipation ratio (FDR) and is a measure of the 
departure from equilibrium. In the equilibrium regime, 
X is equal to unity and the usual FDT is recovered. In 
the out of equilibrium regime, X is generally less than 
unity. For large values of t and tw , it has been hypothe- 
sized that X becomes a function of time only through 
the correlation function C and a Quasi-FDT (QFDT) 
relation is obtained, 
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An important property of this function X{C) is that 
it provides indirect information on the structure of the 
phase space. The ratio Teff{C) = has been inter- 

preted as an effective temperature which is generally 
larger than the heat bath temperature T. This effective 
temperature has been used to classify glassy systems into 
three main groups: in models with a single pure state 
such as a pure ferromagnet and random Ising systems 
[loj , Tef f (C) is infinite since X = ; in models of struc- 
tural glasses 0, [3 , Lennard- Jones glasses jl^ Q and 
p-spin models in d = 3 , Tg^ f (C) > T is finite; in sys- 
tems with full replica symmetry breaking as in several 
finite dimensional spin glasse s ITo . \v\. and mean field 
(MF) models for spin glasses [l9l, Teff{C) is a nontrivial 
function of C. 

Geometrically frustrated antiferromagnets can also ex- 
hibit glassy behaviour and other novel kinds of low tem- 
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perature magnetic states mil 113 which are quite dif- 
ferent from those observed in conventional magnets. Re- 
cently we have performed a numerical study of the two 
dimensional eas y-ax is Heisenberg antiferromagnet on the 
Kagome lattice |23| by computing its static and dynamic 
properties. The magnetization indicates a finite Tc but 
Monte Carlo snapshots of the individual spins below this 
temperature do not indicate any long ranged spatial or- 
der. Rather, the individual spins appear to be in a frozen 
state similar to a glass. The three spins on each ele- 
mentary triangle form a distorted 120" planar state with 
a net magnetization in the z-direction but there is no 
sublattice order. We extracted the critical temperature 
and the critical exponents associated with the magneti- 
zation, the susceptibility and the correlation length. We 
also studied the two-time spin-spin correlation function 
at high and low temperatures. We have found that non- 
exponential relaxation sets in at a temperature T* > Tc- 
The relaxation time increases according to a power law 
and diverges at a temperature Tg < Tc where a tran- 
sition to a glassy phase is located. Below Tg, we have 
found clear evidence for the presence of aging effects in 
the autocorrelation function from off-equilibrium dynam- 
ics. The aging effects obey the same scaling laws that are 
observed in spin glasses and polymers. 

In this letter we report the first observation of the vio- 
lation of the FDT in this low dimensional geometrically 
frustrated system without disorder. We determine both 
the FDR and the effective temperature Te//(C). Our 
results for X{C) indicate that the system behaves like 
a generalized mean field model of glasses in which there 
is a one step replica symmetry breaking. In the aging 
regime, X{C) has an approximate linear dependence on 
the temperature. 

The Hamiltonian describing the model is, 

H = jY^iS^S^ + SfS] + AS^S]) (3) 

i<j 

where {Sf^a — x,y,z) represents a classical three com- 
ponent spin of unit magnitude located at each site i of 
a Kagome lattice. The number of sites N = where 
L is the linear size of the lattice and the exchange inter- 
actions arc restricted to nearest- neighbour pairs of sites. 
The parameter A describes the strength of the exchange 
anisotropy. In the following we restrict our consideration 
to the value A = 2 which corresponds to the easy axis 
case. 

In order to study the non-equilibrium properties of this 
easy axis frustrated system, we calculate the spin-spin 
autocorrelation function [2^ 

C{t + t^,t^) = < S^it + t^)S^^it^) > (4) 

i 

where < ... > means an average over thermal histories 
and tw is the waiting time measured from some quench- 



ing time to = 0. A second quantity of interest is the as- 
sociated response function to an external magnetic field 
hi(t) applied at 

M{t + t^,t^) = J R{t + t^,s)h{s)ds. (5) 

The off equilibrium susceptibility x{t + tnniw) can ob- 
tained by dividing by the field magnitude and the QFDT 
relation (2) allows us to write for long times (i, — > 00) 

Txit + t^,tw)^^M{t + t^,t^,)= [ " X{C)dC. 

(6) 

If the FDT is satisfied {X{C) = 1), then we obtain 
the linear relation Tx{t + tyj,tyj) — Cq — C{t + t^,t^). 
A departure from this straight line in a Tx vs C para- 
metric plot indicates a violation of the FDT and yields 
information about the function X{C). In Ising systems, 
the upper limit Co in (6) is unity but in this Heisen- 
berg system the value depends on both T and the value 
of the easy axis anisotropy A. At T = it is given by 
^0 ^ ''^t(A+i)'^ must be determined numerically at 
finite T. 

For each Monte Carlo run, the system is initialized in a 
random initial configuration corresponding to a quench- 
ing from infinite temperature to the temperature T. We 
then allow the system to evolve for a time and we 
then make a second copy of the system. Using the orig- 
inal copy we compute C{t + tw,tw) as a function of the 
observation time < > 0, for different values of and T. 
In the second copy we apply a perturbation in the form 
of random small magnetic field h{i) = /i^e,; in order to 
avoid favoring one of the different phases The 
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FIG. 1: Correlation function, C{t + tui,tui) ( upper curves), 
and the integrated response function, Tx{t + t^jt^), (lower 
curves) at T = 0.04 and /i^ = 0.01 for waiting times = 
100, 5000 
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FIG. 2; The integrated response vs the correlation function at 
T = 0.04 for waiting times t„ = 100, 500, 1000, 5000 and field 
strength = 0.01 with L = 60. The straight line corresponds 
to Tx = 0.617 - C. 

are taken from a bimodal distribution (e; = ±1) and the 
strength of the field ft.^ is taken to be small (typically 
between /i^=0.002 to 0.01) to ensure linear response. We 
compute the staggered magnetization in this second copy 

M{t + t^,U,) = ^ ^ < Sfit + t^)e, >, (7) 

■i 

whose conjugate field has magnitude and where the 
over line means an average over the random variables e^. 

Typical data for the integrated response and the corre- 
lation function are shown in Fig. 1 where both functions 
exhibit a stationary part for at short times t << and 
an aging part for t >> tw In Fig. 2 , a parametric plot of 
the integrated response versus autocorrelation function is 
shown for T — 0.04 obtained with a perturbing field mag- 
nitude /i^ = 0.01, for waiting times i„ = 100,500,1000 
and 5000 Monte Carlo steps (mcs) . It is clearly seen that 
for short times {t « tu,) TTI and FDT hold and the 
straight line has a slope of —1. For larger times {t » tu,) 
a departure from the FDT line is observed indicating that 
the system has fallen out of equilibrium. A closer look 
at the data in the aging regime shows that , while the 
spin-spin autocorrelation function, C{t + t^jt^) is de- 
caying to zero as i — > oo, the associated response func- 
tion M(t + tw,tu]) keeps growing for all waiting times 
t-u] and its value depends on tyj. These features are re- 
flected in a nontrivial function X{C) . Similar curves 
have been obtained in previous studies of several glassy 
systems 0, 0, Q, such as p-spin models and struc- 
tural glasses. 

The curves in Fig. 2 show a strong dependence on t^,, 
indicating that we are far from the asymptotic regime. 
This dependence on tw could possibly be explained by the 
fact that we are using waiting times that are too small 
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FIG. 3: The integrated response vs the correlation function 
at T = 0.055 for various sizes L, waiting times and field 
strengths . The value of Co = .614 and the two solid lines 
interesect at C = qsA = -55 with slopes —1 and —0.78 . 

[l^ . Indeed, by performing simulations at a slightly 
higher temperature T = 0.055 we see in Fig. 3 that the 
curves clearly indicate: (i) two well separated time scales, 
(ii) no dependence of the FDR on in the aging regime, 
and (iii) no finite size effects. The presence of two time 
scales indicates that the system falls out of equilibrium 
at the point where the two straight lines intersect. If we 
identify this point with the i^, — » oo limit of the two- 
time autocorrelation function, then the FDR X{C) can 
be described 0| in the two regions as follows 

^^^^^ " { X[<1EA -C] + [Co - qEA], C < ITa 

where qea is the value of C at the crossing point and 
plays a role similar to that of the Edwards- Anderson or- 
der parameter in the mean field theory of spin glasses 
[i^ l- Fitting the data in the out of equilibrium regime 
to the straight line of eqn(8), will give us both the FDR 
X{C) and qea- 

The values of X{C) and qea obtained from this pro- 
cedure are displayed in Fig. 4 as a function of T us- 
ing the data for = 2000 and L = 60. We can see 
that at high temperatures X is equal to 1 and the usual 
FDT holds, whereas for T < T^ff our results are well 
approximated by the straight line X{T) — T/T^ff with 
Te// = 0.074 ± 0.001 . As the temperature increases the 
intersection point qEA decreases to zero with a sharp drop 
at Teff for this lattice size L = 60. Our previous results 
[isj for this model indicated a glass transition at a slightly 
lower temperature Tg — 0.071 ± 0.002 where the relax- 
ation time diverged. One would expect that T^ff > Tg 
for systems with finite values of L and and that it will 
approach Tg for larger system sizes and longer waiting 
times. An extension of the present work to larger val- 
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FIG. 4: The quantities X and qea for = 2000 as functions 
of the temperature. The straight line is the prediction of the 
approximation X{T) = T/Teff. 



ues of both L and tw will also help to determine whether 
QEA approaches zero continuously or discontinuously at 
the glass transition. 

In summary, we have shown that the fluctuation- 
dissipation ratio can be computed in this disorder free 
Kagome easy axis antifcrromagnet. Several nontrivial 
features predicted by the mean field theory spin glasses 
are also present in this model. Namely, long-ranged mem- 
ory effects and a constant FDR X{C). These features 
correspond to a phase space structure similar to that 
of spin systems undergoing a one step replica symme- 
try breaking. The temperature dependence of the FDR 
X{T) varies linearly with temperature X{T) = T i Tpjf 
as observed in binary mixtures of soft spheres |l3| and 
Lennard- Jones glasses ■25]. All of our results are con- 
sistent with the existence of a glassy state at low tem- 
peratures in this geometrically frustrated model. The 
ideas developed for mean field descriptions of glassy sys- 
tems seem to have an application far beyond the original 
models. Further experimental studies of these low di- 
mensional frustrated magnets could provide some insight 
into the behaviour of glassy materials. 
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